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I. INTRODUCTION 



Recent work has seen unprecedented advances in the preparation of Bose-Einstein condensates (BEC) of dilute 
alkah vapors . These quantum degenerate systems have paved the way for numerous innovative studies of weakly 
interacting Bose gases. An important area within this field is the generation and study of quantized vortices on 
atomic mesoscopic rings and the potential to study effects due to persistent currents and Josephson effects p|-p^. 
Central to such studies is investigation of geometries for generation of a BEC in a toroidal trap ||ll|,|l2|. Typically 
BEC is created in a magnetic trap that confines atoms in weak-field seeking states and is thus dependent on the 
atomic hyperfine states. A toroidal trap can then be formed by piercing a magnetic trap with a blue-detuned laser 
at its center [Q. However, the state dependence of the trapping is a limitation for advanced studies including multi- 
component spinor condensates. Further, the magnetic field is primarily dictated by the trapping requirements setting 
limitations on its s pat ial form and amplitude. Optical dipole traps can potentially circumvent these problems being 
state independent . Importantly, the spatial form of the optical trap is dictated by the light beams used. This 
allows one to potentially generate arbitrary shapes of condensate p^ . Experimental work has shown the ability to 
transfer a condensate from a magnetic trap to an optical dipole trap created from a tightly focused Gaussian light 
beam @. 

The circularly symmetric Laguerre-Gaussian (LG) laser modes have generated substantial interest in recent years. 
This stems from the identification that they possess an orbital angular momentum of £h per photon . This is in 
addition to the spin angular momentum associated with the polarization state of the beam. The azimuthal index £ 
refers to the number of 27r phase cycles around the mode circumference. A given mode has p -I- 1 radial nodes in the 
mode profile where p is the other index used in the LG mode description. One finds that LG beams with radial index 
p — are in the form of an annulus | f6 1 . There are several techniques for production of LG modes including the use 
of a cylindrical lens mode converter [17| and holographic methods |l8|. The latter technique allows generation of LG 
modes from the fundamental output of a laser beam. The annulus becomes thinner as one increases the azimuthal 
index £ of the mode [T^] . 

In this work we present a technique for generating a toroidal BEC employing a Laguerre-Gaussian light beam. 
Single-ringed forms of such beams offer for the first time a direct possibility of generating ring-shaped condensates in 
an optical potential. We study the operation of this dipole trap as a function of azimuthal mode index and discuss 
the loading of condensates into such optical traps from a variety of initial conditions. 



II. BASIC MODEL AND EQUATIONS 

A. Gross-Pitaevskii equation 

The basic model we consider is shown in Fig. 1 and comprises a two-dimensional (2D) BEC whose normalized 
mode profile is frozen by tight confinement along the z-direction, and which is placed at the focus of an off-resonant 
LG laser beam of frequency ujl incident along the z-axis that provides a two-dimensional optical potential ]l3[ | . The 
tight confinement may be provided, for example, by an independent sheet potential using a scanned optical dipole 
potential |T^. Then the Gross-Pitaevskii equation (GPE) describing the macroscopic wavefunction ^(r^,t) for the 
quasi-2D motion can be written as ||f9| , |20|] 

tn^ = -—Vlij + U{r^)i^ + gN\^\^^, (1) 

where rj^ is the two-dimensional position vector in the plane perpendicular to z, the corresponding two-dimensional 
Laplacian, M is the atomic mass, and g is the effective short-range interaction strength. The potential term on the 
right-hand-side 

describes the 2D optical dipole potential with A = cj^ — the laser detuning from the optical transition frequency w^i, 
r the natural linewidth of the optical transition, /gat is the resonant saturation intensity, and /(rj_) = ^eQcn\E{r±, 0)^, 
assuming the BEC is centered at z = and has a thickness less than the Rayleigh range of the focused beam. In 
addition, we assume that the intensity profile of the LG beam is undistorted upon propagation through the 2D BEC, 
which is reasonable for a thin BEC. We remark, however, that the LG beam becomes spatially phase-modulated upon 



2 



transiting the BEC, resulting in a changed far field profile of the beam that could prove a useful diagnostic of the 
BEC density profile. 

For later purposes we point out the GPE conserves both the wave function norm 
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and also the effective single-particle Hamiltonian 

H{t) = J (fr 
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B. Toroidal optical dipole potential 



Our goal in this work is to explore theoretically the use of LG beams for producing toroidal optical dipole traps, the 
dipole potential being proportional to the laser intensity in the limit of large detunings considered here. The intensity 
profile of an LG beam at its focus {z = 0) in cylindrical coordinates (r^ = (r, 9)) takes the form 
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where Pq is the beam power, £ is the azimuthal mode index, the field having a variation exjp(i£9), p is the radial mode 
index which is the number of radial intensity maxima, Wp/ is the mode spot size, and Lp is a generalized Laguerre 
polynomial. Here we create a single-ringed toroidal trap using a red-detuned laser field (A < 0), so the atoms are 
light-seeking, using an LG beam with p = 0, in which case = 1, and we hereafter drop the p index and take £ > 0. 



Then for each value of £ the intensity profile has a single maximum at = | ]21| ]. Thus, in order to produce a 
toroid of a fixed radius ri — vt, we need to choose the spot size Wi for each azimuthal mode such that 
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Furthermore, for a given £ the peak intensity at the toroid center r — rx may be written as 



2Pn 



£\ 



h 



£' 



£\ 



(6) 



(7) 



with Ii = e ^Pq/ttv'^ the peak intensity for £ = 1 |21|. Here in the expression for we have used Stirling's formula 
£\ ^ V2tl£ ■ £^/e^, and we see that the peak intensity scales as VI. 

Gathering the above results together, for a fixed laser power Pq we may write the optical dipole potential (eq. 
in the following form which is useful for comparison between different values of £ 
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Ui being the optical dipole potential well depth for ^ = 1 at the toroid radius tt- Note that Ue is negative for a 
laser red detuned from resonance. Figure 2 shows the normalized optical dipole potential Ui{r)/\Ui\ versus r/rx for 
various azimuthal mode indices £, a red-detuned laser, and fixed laser power and toroid radius, i.e. fixed Ui. Here 
we see that as £ increases the toroidal trap becomes deeper and tighter, meaning that higher-order LG beams present 
advantages for making tight toroidal dipole traps. 
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III. GROUND STATE TOROIDAL SOLUTIONS 



In this section we elucidate the ground state properties of toroidal traps using the harmonic oscillator and Thomas- 
Fermi approximations. These approximate solutions amply illustrate the key features of the problem. 



A. Thomas-Fermi solution 



Setting 'ip{r,9,t) — (j)(r) exp(—ifit/h) for the cylindrically symmetric ground state we obtain 
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(10) 



with fi the chemical potential. To obtain insight into the solutions of this equation we consider the harmonic oscillator 
approximation to the toroidal dipole potential (eq. m)) trap around r — rx, 
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where Ui is negative for a red detuned laser, as assumed throughout the paper. If we furthermore employ the Thomas- 
Fermi approximation [p2|j2^ ], in which the kinetic energy term is neglected in comparison to the mean field energy, 
then we obtain the approximate GPE 
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where the effective harmonic oscillator frequencies are 



where Stirling's formula was used in the last expression. Equation (|lz 
r < rx + Ar^) 
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has the approximate ring solution {rx — Ar^ < 
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where the ring width is given by 
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These solutions correspond to ring-shaped BEG density profiles which are peaked at r = and have a width 2Ar^. 
We see that as the azimuthal mode index £ is increased with all other parameters fixed the radial harmonic oscillator 
frequency increases fie cx £3/4^ and the ring width decreases Are oc £ Thus, using higher-order LG beams provides 
an advantage for making tight bound toroidal rings as commented earlier. 

The 2D Thomas-Fermi solution above should be applicable under the combined conditions gN >> /2M and 
Are/rx << 1, which ensure that the mean-field energy is greater than the kinetic energy, and the ring width is less 
than the toroid radius. The latter condition is also required for the validity of the harmonic oscillator approximation. 



B. Scaling 

To explore the parameter space for the toroidal trap it is useful to introduce appropriately scaled units. For a 
trap of radius rx we scale all lengths as p±^ = t^/tx, all energies are scaled to Tiiox — h^/Mr'^, and we introduce a 
dimensionless time r = ujxt- In these units the GPE becomes 

i^ = -lvl^ + ue{p)^ + TrvN\^\''ip, (16) 
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where (p = rTip, ug = Ui/TujJt, and rj = Mg/irh'^ is a dimensionless measure of the repulsive many-body interactions 
^ , ^J24| . We shall use the above scaled GPE in our numerical simulations of loading toroidal optical dipole traps. 

Turning to the Thomas-Fermi solution, in dimensionless units the radial harmonic oscillator frequency and the ring 
width become 
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and the chemical potential is given by 



A77 



1/3 



1 -f 2 • sgn(ui) 



■qN 



^-1/2 



2/3' 



(17) 



(18) 



where Stirling's formula has been used. These solutions depend on the dimensionless well depth ui for (, — 1 which is 
given explicitly by 



Ui 



'MPn 



huT ShA \ Tr/sat 



(19) 



and is independent of the toroid radius. 

In dimensionless form the condition gN >> /2M for applicability of the Thomas- Fermi approximation becomes 
11N >> l/(27r). Combining this with the condition Ari/rx < 1 yields the constraint on the number of particles 
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which shows that traps using higher azimuthal mode index can hold more atoms, all other parameters being equal. 



C. Parameter values 



To illustrate the basic scales involved in toroidal traps we consider the case of the Xa — 589 nm transition of Na 
using a laser wavelength of Xl = 985 nm and a toroid radius of rx = 10 /im, these values being used throughout this 
paper. Then using the parameter values /sat = 63 W/m^, F = 27r x 9.89 MHz, we obtain ue = uiVI = —253 x PqV^j 
with Po in milli-Watts. Using hujr = 2.1 x 10^^" K, with corresponding time scale 1/ujt = 36 ms, for a given 
azimuthal mode index £ this potential depth translates to ui — —53 x Pq^/IuK. The corresponding effective oscillator 
frequencies in the toroidal trap are then fli = 2Trx 0.14Pg^^^i'^/^kHz. Consider therefore the case of Pq = ImW so that 
Ml = —253. For two-dimensional BECs a characteristic value for the many-body parameter is 77 « 10"'^ Then 
for N = 10^ and £ = 6 we have Arg = 1.9/zm, and fig — 2Trx 0.54kHz. This is not a particularly tight ring but tighter 
confinement can be achieved by increasing the azimuthal mode index. Figure 3 shows the variation of the oscillator 
frequency f2^/27r in kHz (solid line), and An in microns (dashed line) both versus i for ui = —253,riN = 100, and 
displays the advantage of using high azimuthal mode index. We note that it is difficult to produce LG modes with very 
high azimuthal index (£ > 6) due to decreasing mode purity from holographic generation. However, other techniques 
exist for generating narrow annular beams, e.g. using an axicon [Q. Typically the radial spread of such beams is 
very small, thus they would allow one to explore the region of tight radial confinement seen for LG beams with £ > 20 
(see Fig. 3). 



IV. LOADING A TOROIDAL OPTICAL DIPOLE TRAP 



Having established that LG beams can provide toroidal traps for 2D BECs the issue arises of how to load atoms 
into the ground state of the trap. It does not seem feasible to condense directly in the toroidal trap as it does not lend 
itself to evaporative cooling and the laser field could lead to heating. The solution is then to load into the toroidal 
trap from an existing BEG, e.g. a magnetic trap or sheet dipole potential trap |l^. Here we consider loading of a 
toroidal trap from an initial 2D ring BEG, for example using a magnetic trap with a blue-detuned laser piercing its 
center and also from a centrally peaked BEG as is the case for conventional harmonic traps. The calculations 
presented here demonstrate that it is in principle possible to load toroidal traps and we use illustrative examples 
which highlight the issues involved. 
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A. Types of initial condition 



We have numerically solved the GPE (|l^) for a variety of initial conditions with the exact optical dipole potential 
in Eq. (||) . We consider the situation that for t < the BEG is in the ground state of a prescribed potential which 
is turned off for i > and the toroidal LG trap is turned on. To model this situation we have considered a variety 
of initial conditions at t = and their subsequent evolution. Rather than dwelling on a specific potential model for 
t < we consider the super-Gaussian initial macroscopic wave functions 

V'(r, 0) = AAe^('-''p-><)'"/^"'" , (21) 

where A/" is a normalization constant, w is the width of the initial wave function, m > 2 is the order of the super- 
Gaussian, m = 2 being the usual Gaussian, and Tpoak is the displacement of the density peak away from the origin: 
for Tpeak = we have a centrally peaked initial density, Tpcak = tt gives an initial ring BEG with its peak at the 
toroid radius. As m increases the initial condition becomes more top-hat like, representative of the broadening due 
to repulsive many-body effects. For the numerics presented here we set m = 8. 

We have numerically solved the GPE (|l^) for a variety of initial conditions with the exact optical dipole potential 
in Eq. (^) subject to the above initial conditions using the split-step beam propagation method |2^]. In implementing 
this scheme we have included absorbing boundary conditions at the edge of the numerical grid, thereby simulating 
losses due to radially outward going atoms. We tested that the numerical results presented were not sensitive to the 
placement of the absorber. Thus, in our simulations the norm n{t) in Eq. (^) of the macroscopic wave function is 
not conserved at unity. However, the norm n{t) < 1 is actually a measure of the fraction of initial atoms that are 
captured in the toroidal trap. 

For most of the simulations presented we chose £ — 6 as this is characteristic of the LG beams that can be generated 
reliably experimentally at present |p^ . 

B. Initial ring BEC 

First we will discuss the loading of the toroidal optical dipole trap with a BEG that is already ring shaped. Although 
the BEC has already the desired shape, this transfer from a magnetic-based ring trap to a toroidal trap is still of 
importance as this trap can confine multi-component condensates in different Zeeman sublevels. For an initial ring 
BEG the ground state of the toroidal trap can be excited by reasonably matching the initial and ground state wave 
functions. Figure ^ shows two examples of the computed dynamics for an initial ring BEG formed in Na with 
Ui = —25.3 {Pq = 0.1 mW), and rjN = 50. For these parameters and i = 6 the Thomas- Fermi theory of the last 
section gives a ring width of Are = 3.2 /j,m. In Fig. ^(a), which shows a gray-scale plot of the 2D atomic density 
\(p{x, 0, t)p versus time in ms along the horizontal axis and radius x in microns along the vertical axis, we used £ = 6, 
1/7 = 4 fj,in, and after a transient the density profile settles down to a steady ring of smaller width than the input. 
The transient involves expansion of the atoms outward from the ring and also inwards towards the origin, and the 
peak at the origin seen for i « 3.2 ms is due to interference between the inwardly propagating circular atomic waves. 
The outward propagating atomic waves leads to an effective loss mechanism which allows for the damping of the 
initial wave function towards the ground state toroidal BEG. The final ring contains n « 90% of the initial atoms. 
Figure |^(b) shows a gray-scale plot of the 2D atomic density for the same parameters as^(a) except £ — 2, for which 
Ar2 = 5.6 /im, so that the input with w = 2.5 /im is now half the expected ring width. After a considerably larger 
transient than in Fig. ^a) the density again settles down but now to a larger ring than the input containing about 
n 80% of the initial atoms. Notice, however, that for t > 15 ms the atomic density shows a further transient, 
and this is caused by weak atomic waves caught around the origin that start to leak out, as revealed by detailed 
examination of the data. 

For these simulations we have used a rather low laser power Pq = 0.1 mW, so that the toroidal trap depth is only 
about 13 nK for £ = 6. We have chosen such a loosely bound toroidal trap as an illustrative example of loading since 
tighter traps (larger |ui|) are generally easier to excite by suitably matched initial ring BEGs as the ground state 
and next excited state have a larger energy separation. Thus the present example shows that there is considerable 
robustness in loading from initial ring BEGs, as expected intuitively. One might expect such shallow traps to be very 
susceptible to heating and losses induced by noise and background gas collisions. However, the heating due to position 
and intensity noise increase with the fourth and second power of the trap frequency fig, respectively, and heating due 
to background gas collisions is proportional to the square of the trap depth |^,^ . Therefore our shallow traps (with 
low trap oscillation frequencies and small trap depth) should be relatively immune to these detrimental effects. 
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C. Initial centrally peaked BEC 



Figure |^ shows representative 2D atomic density plots from simulations of loading from an initial centrally peaked 
(''peak = 0) Na BEC with m — 8,ui — —100 (Fq = 0.40 mW), the corresponding trap depth being 21.2-\/£nK, and 
rjN = 400. The plots are for different initial BEC widths (a) w = = 10 /im with i = 6, (b) w = lA with 
£ — 6, and (c) w — lA /zm, with £ = 2. Case (a) shows very little signs of trapping in the toroid which is centered 
at rx = 10 fim, but (b) with an increased BEC width of w = 14 /im now shows substantial trapping with N — 90% 
of the initial atoms trapped. What distinguishes these two case is that the initial Hamiltonian in Eq. (^) , which we 
evaluated numerically, is H{0)/hujT — 148 in case (a) and H{Q)/'hu)T = —16.6 for case (b). The Hamiltonian has 
three contributions, the kinetic energy which is positive, the optical dipole potential which for a red-detuned laser is 
negative, and the nonlinear term due to many-body repulsion which is positive. As we simulate the effects of atom 
losses due to outward propagating atoms using the absorbing boundary, the atom number decreases and so typically 
does the Hamiltonian as the initial BEC evolves. Physically, as we consider red-detuned traps that have negative 
optical dipole potentials which go to zero away from the trap, atoms with energies above the trap energy of zero will 
tend to be lost from the toroidal trap. For case (a) the initial Hamiltonian or energy is above the trap energy and the 
atoms fly above the toroidal trap under the influence of the repulsive many-body effects. Indeed the initial density 
proflle with w = 10 /im overlaps the toroidal trap very little, leaving only positive contributions to the Hamiltonian 
in Eq. (Q). In contrast, for the wider initial BEC in Fig. ^(b) the initial density profile overlaps the toroidal trap 
giving a net Hamiltonian less than the trap energy, so trapping becomes possible. We remark that the requirement 
that the initial Hamiltonian be less than the trap energy is a necessary but not sufficient condition for trapping, as 
there is still dependence on the initial density proflle. For example, Fig. ||(c) is the same as in ||(b) except ui = —253 
(Pq = 1 niW), the corresponding trap depth being 130 nK for £ = Q. In this case H{0)/hu!T — —216, well below the 
trap energy, but Fig. ||(c) shows that the atoms do not cleanly load into the ground state of the toroidal trap, but 
rather the density displays undamped oscillations on the time scale of the simulation. This oscillatory behavior is 
typically a feature for deeper traps. We have begun to explore methods of smoothing the loading of deeper toroidal 
traps, such as tapering the turn-on of the optical potential and phase-imprinting the initial BEC so that it moves 
towards a ring. These methods do smooth the loading process and we shall report on these issues along with the 
effects of beam misalignments in a future publication. 



V. SUMMARY AND CONCLUSIONS 



In this paper we have shown that Laguerre-Gaussian beams provide a flexible means for forming toroidal optical 
dipole traps in 2D atomic BECs, and that the toroidal traps can be loaded from initial conditions representative 
of conventional magnetic traps. It remains to be seen how tight the atoms can be confined in these rings, and 
detailed numerical studies are underway. This work is a first step towards developing toroidal traps which shall act 
as mesoscopic rings for atoms for a variety of basic and applied studies P-|l0||. For example, once a BEC is prepared 
in its ground state of the toroidal trap a vortex state of variable angular momentum can be excited using Raman 
coupling (involving a second LG beam) to another Zeeman sublevel [ p9[ , hence allowing for studies of ring vortices and 
persistent currents on a torus iQ,^,^- Also, as our toroidal trap does not involve a magnetic trap it allows for studies 
of multi-component BECs trapped on a ring ||l^. In the limit of small number of atoms there is also the possibility 
of realizing a Tonk's gas of impenetrable bosons on a ring |^-^ , which has been predicted to exhibit dark solitons 
p3[ . Furthermore, by looking at higher-order LG modes with different radial index p, and hence multiple concentric 
rings, we can create coaxial toroidal traps for the BECs, hence allowing for radial tunneling between condensates. 
In particular, multiple rings could create a circular grating for atoms which could in principle act as the feedback 
mechanism for a 2D atom laser, by analogy to circular grating optical lasers On the applied side, once a ring BEC 
is formed we can further pierce it with blue-detuned lasers at positions along its perimeter to form tunnel-junctions 
akin to superconducting links. One can then envision a range of sensors, e.g. rotation, based on the sensitivity of the 
tunneling current to any perturbation of the system. In principle, large rings could be made by first loading a small 
ring and then adiabatically expanding the LG light mode incident on the BEC, hence very sensitive matter- wave 
sensors of inertial forces could result. 
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Army Grant No. DAAD 19-00-1-0169. KD acknowledges the support of the UK Engineering and Physical Sciences 
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FIG. 1. Our basic model comprises a two-dimensional BEG which is illuminated by a red-detuned LG beam traveling along 
the «-axis and which acts as a toroidal trap. 

FIG. 2. Normalized optical dipole potential Ui{r)/\Ui\ versus r/rr for various azimuthal mode indices £ for fixed laser power 
and toroid radius. For the red-detuning assumed here, as £ increases the toroidal trap becomes deeper and tighter. 

FIG. 3. Variation of the oscillator frequency Qi/2tt in kHz (solid line), and width Are in microns (dashed line) both versus 
£ for ui = -253, ryiV = 100. 

FIG. 4. Loading of a toroidal trap from an initial ring BEG with rpeak = 10 /im and m = 8: Gray scale plots of the evolution 
of the atomic density \(p\^ for ui — —25.3, rjN — 50, and (a) w — 4 fim with £ = 6, (b) w = 4 with £ — 2. 

FIG. 5. Loading of a toroidal trap from an initial super-Gaussian with m — S: Gray scale plots of the evolution of the atomic 
density \ipf{x,0,t) for ui = -100, r]N = 400, £ = 6, and (a) w = rr = W with H{fi)/hujT = 148, (b) w = U urn with 
H{Q)/hwT = —16.6. (c) is for the same parameters as (b) except u\ — —253 giving H{Q)/Tiujt ~ —216. 
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Figure 1: E. M. Wright et al., Toroidal Optical Dipole Traps... 
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Figure 2: E. IVI. Wriglit etal., Toroidal Optical Dipole Traps... 
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Figure 3 
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Figure 4(a) 
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Figure 4(b) 
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Figure 5(b) 
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Figure 5(c) 
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